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Abstract 

Using 13.6 fb _1 of continuum two-jet e + e~ — > cc events collected with the 
CLEO detector, we have searched for baryon number correlations at the pri- 
mary quark level. We have measured the likelihood for a Aj charmed baryon 
to be produced in the hemisphere opposite a A c relative to the likelihood 
for a charmed baryon to be produced opposite an anticharmed meson D; 
in all cases, the reconstructed hadrons must have momentum greater than 
2.3 GeV/c. We find that, given a A c (reconstructed in five different decay 
modes), a A+ is observed in the opposite hemisphere (0.72 ± 0.11)% of the 
time (not corrected for efficiency). By contrast, given a D in one hemisphere, 
a A+ is observed in the opposite hemisphere only (0.21 ± 0.02)% of the time. 
Normalized to the total number of either A c or D "tags", it is therefore 
3.52±0.45±0.42 times more likely to find a A+ opposite a A c than a D me- 
son. This enhancement is not observed in the JETSET 7.3 e + e~ —* cc Monte 
Carlo simulation. 
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I. INTRODUCTION 

At high energy, quark fragmentation can be calculated using Perturbative QCD (PQCD). The 
results are often referred to as QCD showers pf. However, these calculations require, for comparison 
with experimental results, some model or experimental input to take into account the eventual 
hadronization of gluons and qq pairs into actual hadrons of known mass. 

Several models have thus been developed, mostly based on a chromodynamic string connecting 
the initial go an d % |2]-||] , where the subscript indicates that these are the primary quarks formed 
by e + e~ annihilation. When stretched because of the initial momentum of go an d go, the string 
breaks by creation of secondary gg pairs, generating two substrings. The process iterates then for 
the new substrings until the physical hadrons are produced. This last stage is the critical one where 
QCD loses its predictive power. 

The model most frequently used to simulate the process e + e~ — > gog — > hadrons is the 
QCD inspired Lund String Model (LSM) ||, implemented in the Jetset Monte Carlo simulation 
package 0]. In order to solve the problem of generating hadrons of known mass, the Jetset imple- 
mentation of LSM, at each step of its iterative process, splits a gg string into a hadron of known 
quantum numbers and mass and a remainder string having the leftover quantum numbers and 
energy-momentum. When, finally, the remainder string has a mass below a suitable limit, this 
string is made into two known hadrons that, together, carry the leftover quantum numbers. 

In this paper we report on correlations between charm baryons in the fragmentation of a cc 
system produced in the reaction e + e~ — > cc at ^/s ~ 10 GeV, where the annihilation occurs in a 
largely low-Q 2 , non-perturbative regime. In contrast to inclusive single-particle production, com- 
pensation of baryon and charm number is a more subtle aspect of quark fragmentation modeling. 
One obvious question is whether baryon compensation occurs locally (e.g., small rapidity difference 
between baryon and antibaryon) or globally (large rapidity difference). Previous studies of AA 
production at ^^=90 GeV found that in events containing both a A and a A, the two particles 
tended to be produced at very similar rapidities ||. 

Consider the case in which a charm baryon is produced in the first step of fragmentation 
(e.g. e + e~ — > cc; c — > A C A); it is possible that both baryon and quantum charm numbers be 
compensated in the opposite hemisphere (e.g. e + e _ — > cc; c — > A C X, c —* A c X). In the limit that 
the c and c quarks fragment independently (presumably true at some sufficiently high energy) , this 
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type of correlation is not expected. The Jetset implementation of the Lund model, because of the 
mechanism outlined above, does not produce such a correlation. 

In the very reasonable approximation of neglecting string splitting by cc tunneling from the 
vacuum, at y/s ~ 10 GeV all observed charmed particles must contain a primary quark. We take 
advantage of this to discriminate between independent vs. correlated fragmentation models. If we 
assume that primary quarks fragment independently, then the number of times that we find a A+ 
baryon opposite a A c antibaryon in an event (i.e., cos6(Af, A c ) <0, denoted "A+|A C "), scaled 

to the total number of observed A c (denoted " ^gJ4r "), should be equal to the number of times 

A c 

that we find a A+ baryon opposite any other anti-charmed hadron H c , scaled to the total number 



of observed anti-charmed hadrons ( Ac ^ Hc ) . In this analysis, we will check the equality of these 

ri c 

ratios; by comparing ratios in this way, we cancel many experimental systematics. Schematically, 
we are comparing the following event topologies: 



c c 

XT «- — A+ 



vs 



c c 

Specifically, we compare the rate of A+|A C production to the rate of A+|D° and A+\D~ 
production. We reconstruct D 0, s and D~'s through the well-measured decay modes D° — ► K + ir~ 
and D~ — > K + n~ir~ , respectively. A+'s are fully reconstructed in the decay modes A^~ — * pK~ir + , 
A+ — > pK®, A+ — > Air + , A+ — > A7r + 7r _ 7r + , A+ — * pKgir + ir~ ^\ and partially reconstructed through 
A+ — > AXj^| S c 's are also studied; these are reconstructed in £+ + — > A+7r + and S° — > A+7r~. 

Under the assumption of independent fragmentation, we expect the relative production ratios 
to satisfy: 



and, 



A+|A e ^ A+l-D 
Ac" ' DO 



A+|A e . A+\D- 



1 



Ac 



D~ 



1 Charge conjugation is implicit. 

2 According to JETSET 7.3 simulations, 95% of the A's opposite an anti-charm tag are charmed 
baryon daughters (predominantly A c ) after imposing a 1 GeV/c minimum momentum requirement 
on the A (pa >1 GeV/c). 
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II. APPARATUS AND EVENT SELECTION 



This analysis was performed using the CLEO II and the upgraded CLEO II. V detectors oper- 
ating at the Cornell Electron Storage Ring (CESR) at center-of-mass energies \fs = 10.52-10.58 
GeV. For 4.5 fb _1 of the data used for this analysis ("CLEO-II" data ||), measurements of charged 
particle momenta were made with three nested coaxial drift chambers consisting of 6, 10, and 51 
layers, respectively. In a subsequent upgrade ( "CLEO-II. V" [[Hj), the inner tracking chamber was 
replaced with a high-precision silicon detector (corresponding to the remaining 9.1 fb" 1 of the data 
used for this analysis). The entire tracking system fills the volume from r=3 cm to r=l m, with 
r the radial coordinate relative to the beam (z) axis. This system is very efficient (e >98%) for 
detecting tracks that have transverse momenta {pt) relative to the beam axis greater than 200 
MeV/c, and that are contained within the good fiducial volume of the drift chamber (| cos 9\ <0.94, 
with 6 defined as the polar angle relative to the beam axis). This system achieves a momentum res- 
olution of (5p/p) 2 = (0.0015p) 2 + (0.005) 2 (p is the momentum, measured in GeV/c). Pulse-height 
measurements in the main drift chamber provide specific ionization resolution of 5.0% (CLEO II. V) 
or 5.5% (CLEO II) for Bhabha events, giving good K/ir separation for tracks with momenta up to 
700 MeV/c and separation of order 2a in the relativistic rise region above 2.5 GeV/c. Outside the 
central tracking chambers are plastic scintillation counters, which are used as a fast element in the 
trigger system and also provide particle identification information from time-of-flight measurements. 

Beyond the time-of-flight system is the electromagnetic calorimeter, consisting of 7800 thallium- 
doped Csl crystals. The central "barrel" region of the calorimeter covers about 75% of the solid 
angle and has an energy resolution which is empirically found to follow: 

a *<%)= 0^ + 1.9-0.1^; (1) 



E is the shower energy in GeV. This parameterization includes noise effects, and translates to an 
energy resolution of about 4% at 100 MeV and 1.2% at 5 GeV. Two end-cap regions of the crystal 
calorimeter extend solid angle coverage to about 95% of 4-7T, although energy resolution is not as 
good as that of the barrel region. The tracking system, time-of-flight counters, and calorimeter are 
all contained within a superconducting coil operated at 1.5 Tesla. An iron flux return interspersed 
with proportional tubes used for muon detection is located immediately outside the coil and in the 
two end-cap regions. 

The event sample used for this measurement is comprised of 9.1 fb _1 of data collected at 
the T(4S) resonance and 4.5 fb _1 of data collected about 60 MeV below the T(4S) resonance. 
Approximately 18 x 10 6 continuum cc events are included in this sample. Charged track candidates 
for protons, kaons or pions must pass the following restrictions: 

(a) The track must have an impact parameter relative to the estimated event vertex less than 
5 mm in a plane perpendicular to the beam axis (f — (ft) and no more than 5 cm along the beam 
axis. The estimated event vertex is obtained by averaging the e + e~ interaction point over a full 
run. 

(b) The track has specific ionization information consistent (at the 99% confidence level) with 
its assumed particle identity. 

(c) The track must have momentum greater than 100 MeV/c. 

All reconstructed charmed hadrons must have momentum greater than 2.3 GeV/c to ensure 
that there is no contamination from £>-meson decays to charm. 
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A. Production Ratios 



We define the single tag yield to be the number of reconstructed events containing one anti- 
charmed hadron H c . This yield (Table |[) is determined by fitting a double-Gaussian signal over 
a smooth, low-order polynomial background functional. The signal function is the sum of two 
Gaussian functions, one narrow and one broad, and is a better representation of the expected 
signal line shape than a single Gaussian function. 

TABLE I. Single Tag Yields obtained from one-dimensional fits (charge conjugate modes are 
implied); statistical errors only are shown. 





Data 


Monte Carlo 




(~60 M hadronic events) 


(~160 M hadronic events) 


A c (5 modes) 


70199 ± 1604 




A c (pK~7r + + pK® only) 


56110 ± 1776 


130380 ± 1943 


£0 _ A+7T- + £++ -» A+7T+ (sum) 


3804 ± 185 


6522 ± 182 


A — » p7T~ 


735343 ± 1198 


2136997 ± 1674 


D° -► K--K+ 


352294 ± 1668 


1045776 ± 2294 


D + -> K-7T+7T+ 


273597 ± 2148 


705357 ± 3232 



The number of double tags (Table |l]) is defined as the number of events in which two specific 
particles are reconstructed in opposite hemispheres (i.e., the opening angle between the two particles 
must be greater than 90 degrees). The total number of double tags for A+|A C , A+|A, A+|D° and 
A\D° is extracted from two-dimensional invariant mass plots, shown in Figures |], ||, and |] 
respectively. 

The total correlated double-tag yield is first determined by performing a sideband subtraction in 
the two-dimensional invariant mass plot. Consider, for example Fig. [l]. We take one-dimensional 
projections of three slices in the candidate A c recoil invariant mass "M reco ^" - the A c signal 
region: (\M recoi i — 2.286| <0.03 GeV/c 2 ) and the two A c sideband regions: (0.04 < \M recoi i — 
2.286| <0.07 GeV/c 2 ). We then subtract the A+ distributions from the A c sidebands from that 
of the A c signal region. Figure |E|. shows the signal distribution and the sum of the two sideband 
distributions. In performing these fits, the double-Gaussian signal shapes are constrained using the 
parameters determined from fits to the single-tag sample. 

We also perform a single fit in two dimensions to extract the signal yields. In this latter fit, a 
two-dimensional Gaussian signal is used to parametrize the peak region, two single Gaussians are 
used to fit the "ridges" away from the peak region (corresponding to true charmed hadrons along 
one axis in association with combinatoric background on the other axis) and a two-dimensional, 
smooth polynomial is used to parametrize the background. Differences in i) the selection of the 
sideband /signal regions and the regions projected, as well as ii) the difference between the sideband- 
subtracted yield compared to the two-dimensional, single-fit yield give a measure of the fitting 
systematics for each measurement. 

We define the "production rate" (Table ID) as the percentage of times we find one specific 
particle in an event opposite a given tag (i.e., the number of double tags, given in Table H divided 
by the total number of single tags, presented in Table 0). Production ratios are formed by taking 
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3031200-037 




FIG. 1. Double tag invariant mass plot of A+ candidates plotted vs. invariant mass of A c 
candidates (A+|A C ) from data. Shown is the sum of the modes: A+ — > pK~ir + , A+ — > pK®, 
A+ — ► A7r + , A+ — > A7r + 7r _ 7r + , and A+ — > pKgiT + 7r~ (and their charge conjugates, in the case of 
A c reconstruction). Horizontal and vertical solid lines designate signal bands. 
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FIG. 2. Double tag plot of A|A C (plus charge conjugate) from data. The A c is selected as in 
the previous figure; the A is reconstructed in A — > pir. Horizontal and vertical solid lines designate 
signal bands. 
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FIG. 3. Double tag plot of A+|D° (plus charge conjugate) for data. Horizontal and vertical 
solid lines designate signal bands. 
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FIG. 4. Double tag plot of A|D° (plus charge conjugate) for data. Horizontal and vertical solid 
lines designate signal bands. 
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FIG. 5. Projections of Figure 1 onto the candidate A+ mass, after requiring that the recoil mass 
be consistent with the A c signal (solid histogram, with dotted curve fit overlaid) and sideband 
(dashed) regions (see text for details). 



TABLE II. Double Tag Yields obtained from two-dimensional fits (statistical errors only). 





Data 


Monte Carlo 


A+|A C 


253 ± 37 


97 ± 40 


At\D° 


722 ± 55 


1905 ± 78 


&t\ D ~ 


556 ± 71 


1281 ± 96 




34 ± 14 


13 ± 9 


A+|A 


1355 ± 72 


1079 ± 85 




29 ± 13 


88 ± 13 


£ c |£>- 


32 ± 18 


59 ± 14 


S C |A 


122 ± 19 


49 ± 12 


A\lf 


2400 ± 69 


8547 ± 117 


A\D- 


2132 ± 90 


6397 ± 144 
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ratios of production rates (Table |TV|) . Note that, to compare to — =^ , we must take into 

account the fact that the latter ratio implicitly includes both sign combinations (— =— plus — ), 

while there is only one unique way to form Ac _\^ . To compare these two rates and form production 

Ac- 



ratios (Table IV), the production rates for A c -tags presented in the table have therefore been 



multiplied by this factor of two. 

TABLE III. Production Rates; statistical errors only. 



Double tags 
Single tags 



Data Fraction Monte Carlo Fraction 



A, 



2 x A+jA c ^ 1Q ± im j x 10 _ 3 ^ 49 ± Q 62 ) x 1(T 3 



-3 



(2.05 ± 0.16) x 10" 3 (1.82 ± 0.08 ) x 10 

(2.03 ± 0.26) x 10~ 3 (1.82 ± 0.14 ) x 10~ 3 

(19.3 ± 1.1) x 10~ 3 (8.28 ± 0.66) x 10~ 3 



£^£_ (0.49 ± 0.20 ) x 10~ 3 (0.102 ± 0.066) x 10" 3 



(0.082 ± 0.036 ) x 10~ 3 (0.068 ± 0.008) x 10~ 3 

(0.067 ± 0.049 ) x 10~ 3 (0.084 ± 0.020) x 10" 3 
^& (32.1 ± 5.2) x 10~ 3 (7.6 ± 1.8) x 10~ 3 

^ (6.81 ± 0.20 ) x 10~ 3 (8.17 ± 0.11 ) x 10~ 3 
(7.79 ± 0.33 ) x 10~ 3 (9.07 ± 0.21 ) x 10~ 3 



TABLE IV. Production Ratios 



Data (stat. and sys. errors) 



Monte Carlo (stat. error only) 



2xA+|A c 


. A+|£>° 


A c " 


' D° 


2xA+|A c 


. A+|D- 


A c " 


D~ 


A|A C 


. M~d° 


a; 


' D° 


A|A C 


. A|D- 


A c " 


• D- 


A|E C 


. a\d° 




' DO 


A|S C 


. A\D~ 


E, 


■ D~ 



(3.51 ± 0.59 ± 0.42) 


(0.82 ± 0.34) 


(3.54 ± 0.70 ± 0.43) 


(0.80 ± 0.35) 


(2.83 ± 0.18 ± 0.26) 


(1.01 ± 0.08) 


(2.48 ± 0.18 ± 0.23) 


(0.91 ± 0.08) 


(4.71 ± 0.77 ± 0.56) 


(0.93 ± 0.22) 


(4.12 ± 0.68 ± 0.52) 


(0.84 ± 0.22) 



From Table III, we see a clear enhancement in the likelihood of producing a A^!~ opposite a A c 
compared to a D. This is observed in the case where the A+ is fully reconstructed, as well as tagged 
inclusively by a A. Note that the fractional enhancement is smaller in the case when the A+ is 
tagged inclusively. This is qualitatively consistent with the expectation that events containing two 
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charmed mesons and two baryons will produce a A\D correlation which will inflate the denominator 

when we construct the ratio: / ^=£. Non-cancelling contributions to the numerator, such as 

_ _ Ac 
AKND\A C , plus decays of charmed baryons other than A c into A, may also be present. 



III. CROSS CHECKS 



A number of cross checks were conducted in order to verify the accuracy of the double tag 
signal extraction. These include signal extractions of two-dimensional "wrong sign" plots (expected 
to have zero signal yield), as well as consistency checks with Monte Carlo and studies of D\D 
correlations. 



A. Null and Wrong sign signals 

Neglecting the doubly Cabibbo suppressed decay D° — > K~ir + , we expect a null signal yield 
from a double tag plot of (mx"^ \ m K-iT+)- Similarly, we expect zero signal yield from a double tag 
plot of (m^-- 7r +m^-+ 7I — |) in the case where the D° and D candidate are in the same hemisphere. 
All of these correlations give signal yields consistent with zero as expected (Table [V]). 

TABLE V. Null Checks (statistical errors only). 



Correlation Yield 

D°\D° (opposite hemisphere) -16 ± 47 

L' ± |L' ± (opposite hemisphere) -10 ± 83 
D°\D + (opposite hemisphere) 55 ± 91 

(same hemisphere) -15 ± 25 

D^D^l (same hemisphere) -116 ± 157 



B. Monte Carlo studies 



The uncertainty in A c production characteristics is expected to be somewhat larger than the 
uncertainty in D° production. We therefore expect the Monte Carlo and data to agree on, e.g., the 
number of D + \D~ double tags per D~ ( — J g ), and the number of D°\D° double tags (Figure |||) 

per D° ( D -^o ), as shown in Table |Vl|.P| 

Within errors, the agreement is good. We have also further ensured that there is no bias in the 
signal extraction due to possible peaking of the background in the signal region, by subtracting all 
the true signal from a Monte Carlo double tag plot and verifying that the measured yield, after 
subtraction of the true generated particles, is indeed consistent with zero. 



3 The measured value of 0.0119, for instance, is qualitatively consistent with the expectation that, 
per D tag, 50% of the time the charm quark will produce a D , which is reconstructed in the 
D° — > K~ir + mode (B(D° — > K~tt + ) ~0.04) with approximately 60% efficiency. 
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FIG. 6. Double tag plot of D°\D°. Horizontal and vertical solid lines designate signal bands. 



TABLE VI. Consistency Check using D\D correlations (statistical errors only). 





Data Raw Rate 


Monte Carlo Raw Rate 


D°\D° yield 


2099 ± 56 


6238 ± 97 


^=0- production rate 


0.0119 ± 0.0003 


0.0119 ± 0.0002 


D + \D~ yield 
D p D production rate 


1260 ± 88 
0.0092 ± 0.0006 


3279 ± 130 
0.0093 ± 0.0004 
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C. Check of expected uncorrelated yield 



Having obtained values for production rates, we can compare these with the values expected 
under the assumption of two hemispheres fragmentating independently. In this (hypothetical) limit, 
the ratio of double-tags to single-tags should be the same as the ratio of single-tags to total charm 
quarks, since each simply expresses the probability of a charm quark to evolve into a particular 
final-state particle. Then, using B{A C — ► X) to designate the branching fraction for A c to decay 
into the final state X, C as the integrated luminosity of our data set, and /( c _>a c ) as the fraction 
of times a charm quark materializes as a A c baryon, we can relate the expected number of single 
tags and double tags as: 

N(Ac),Single Tags = £ " <j{e + e~ -> Cc) ■ /( C -A c ) " #( A c ~> X) ■ 6 Ac ^ X (2) 

and, for the probability that, given a reconstructed A c , we reconstruct an opposite A c 
7V(A+|A C -) 

Double Tags — 

/( C ->A C ) ' ^(^c X) ■ e Ac _^ x ■ N(A c ) S ingle Tags (3) 



In equation (2), for the number of single A+ tags, we use half the total number of observed 
A+ + A c (70199/2). In equation (3), the term /| c _ >Ac ^ indicates that the fraction of times that 
a charm quark produces a A c may be different if the event already contains a A c produced from 
the corresponding anticharm quark. This is, of course, exactly the correlation factor we wish to 
ultimately measure. The prime on e Ac ^x indicates that the efficiency for reconstructing a A c may 
be higher for events in which a A c has already been reconstructed due to geometrical correlations 
- since the charm and anticharm quarks are back to back in an e + e~ — > cc event, reconstruction of 
one charmed particle ensures that the corresponding antiparticle is in a good acceptance region of 
the detector. Dividing equation (3) by equation (2) and multiplying each side by N(A c )singie Tags 
we obtain, for the number of double tags: 

. r , A + nr -s (N(A c ) S ingleTags) 2 ■ /( C -A c ) " £ A C ^X , A . 

N{At\A c ) Doub leTags = ~ r 7 T ~ =x~ T ~ (4) 

C ■ a(e+e -► cc) • /( c _»a c ) ■ ^a c ^x 

The difference between /( c _ >Ac ) and /| C _ >A ) represents the enhancement in production of the A c 

e ' C /' 

when a A c is present. With Ac ^* = e geome tr y and = /correlated (=1 if independent frag- 

mentation holds), we can estimate our expected double tag to single tag ratio as follows: From 

the luminosity, we calculate the number of e + e~ — > cc events using as inputs the e + e~ — > qq 

cross-section (3.3 nb), and taking ^ + ^_~_^£ =0.4. Assuming that /correlated is equal to unity for D\D 

events, we can solve Eqn. (4) for egeometry by comparing the total number of D°\D° and D + \D~ 

double tags we would expect to find (assuming independent fragmentation) with the actual number 

of D°\D° and D + \D~ double tags. This is calculated to be 

(2099)(0.4)(3.3)(13.6 x 10 6 ) lrtrtJ , . . , 

egeometry = (352294/2) 2 ^ {statistical errors only), 

using the total luminosity £=13.6 fb -1 , and the single tag and double tag values for our D° and 
D°\D samples, respectively. For the D + and D + \D~ samples, We obtain egeometry — 

1.21 ±0.09. 
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Using a value of e geometry = 1.21, a similar calculation can be performed for the expected number 
of A+|A C double tags, based on the total number of cc events and the single tag yield, giving an 
expected value of 83.0±2.6 (statistical error only) double tag events. This number is then compared 
to the measured number of A+|A C double tags to obtain an estimate of f correlated = ffg^fi = 
3.04 ± 0.45. Despite the roughness of this approach, the observed enhancement is consistent with 
our calculated production ratios. 



IV. STUDY OF THE A+|A C SIGNAL CHARACTERISTICS 



We have compared some of the production characteristics of the observed A^~|A C signal with 
the A+|D° and A+ \D~ signals. No obvious difference is found in either the A+ momentum spectrum 
(or A c , in the charge conjugate case) or the polar angle distribution of the A+ (see Figs. ^ and|| 
respectively) for the A c -tagged or the .D-tagged samples. We also find that the observed A+|A C 
cross-section is statistically the same for the data taken on the T(4S) resonance as data taken 
just below the T(4S). As expected (Figure |9|) if the charmed baryons are following the direction 
of the original charm/anticharm quarks, the observed opening angle between the A+ and A c 
peaks at 180°. Finally, to investigate the possibility that the A+|A C signal was associated with 
the production of a 4-baryon system, we plot the double tag yield in cases where there are well- 
identified protons found in the same event (Fig. 10). No signal is observed in such a case. 




2.20 



2.62 



3.04 



3.46 



3.88 



4.30 



4.72 



Momentum of A c (GeV/c) 



FIG. 7. A c (and charge conjugate) momentum distribution for the three double tag signals 
A+|A C , A+|.D , and A+|D~. Squares represent the A+|A C signal, upright triangles represent the 
Af\D° signal, and inverted triangles represent the A+|L>~ signal. These three distributions are 
qualitatively similar, indicating similar production dynamics. All distributions have been normal- 
ized to have unit area. 



14 



0.50 r 



■ 0.25 - 



a a; I D° 

v a; i^r 

♦ D°|D° 



< 



p,/p 



FIG. 8. Polar angle distribution with respect to the e + e~ beam axis (p z /p) for the double tag 
signals A+|A C (squares), A+|L>° (triangles), A+|D~ (inverted triangles), and D°\D (diamonds). 
All distributions have been normalized to have unit area. 




0.50 



FIG. 9. Production angle between the A+ and A c (squares) compared with A^\D (upright 
triangle), Af\D~ (inverted triangle), and D°\D° (diamond) event samples. As expected, A+ are 
produced back-to-back with respect to A c tags. All distributions have been normalized to have 
unit area. 
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FIG. 10. Double-tag A+|A C signal in events containing an identified antiproton or proton. 

V. £ c (|)+ PRODUCTION 

We have also investigated the possibility of correlated S C |S C production. Of particular interest 
here is the possibility of observing an enhanced production rate relative to our observed A C |A C 
correlation. Such an enhancement may be indicative of spin correlations at the first rank in frag- 
mentation. S c 's are reconstructed through their decay mode into A c and a soft pion tt s : S c — > A c 7r s . 
For the highest efficiency and best resolution, only charged pions irf- are used in reconstructing E c 's. 
Although statistically limited, a direct search for S C |S C double tags yields 6 events in the signal 



region (Figure 11); the extrapolated background under the signal is 1.8 ± 0.3 events. 

Using the total number of S c single tags observed (3804 ± 185), the total number of A c single 
tags observed (70199 ± 1604), and the total number of A C |A C double tags observed (252.5 ± 37.4), 
we can estimate the number of expected S C |S C double tags (assuming that the S C |S C correlated 
production rate is the same as the correlated A C |A C production rate) as: 252.5x( 7 3 8 1 ° 9 q ) 2 =0.74 zfc 
0.11. We therefore see a statistically very limited indication of S C |S C correlations, suggesting a spin 
correlation at the first step of fragmentation. More data are necessary to elucidate this situation. 



VI. SYSTEMATICS 



We expect most of the systematic uncertainties associated with measurement of both numerator 
and denominator of the production rates to cancel in calculation of the production ratios. The sen- 
sitivity to the definition of "signal" and "sideband" mass regions in the one- and two-dimensional 
subtraction was determined by calculating the yields as we vary these parameters. We assign a sys- 
tematic error of 6% due to "signal parameterization uncertainty" based on this study. We have also 
performed a second (complete) set of fits using two-dimensional Gaussian parameterizations of the 
signal and two-dimensional polynomial parameterizations of the background. The two techniques 
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FIG. 11. Double-tag S C |S C mass difference plot. Horizontal and vertical solid lines designate 
signal bands. 



show excellent agreement, and are typically within 5% of each other, for each two-dimensional 
extraction. When we vary the selection criteria (e.g., particle identification criteria, and minimum 
momentum requirements) used to define our hadrons (A c , D°, and D + ), we observe a maximum 
7% variation in the production ratios we calculate. As an example of our investigation of the mo- 
mentum dependence of the A^~A C production ratio enhancement, Table VII shows the production 
ratio in the two-dimensional (pa c ,Pjj ) space. 

TABLE VII. Production ratio enhancement tabulated as a function of A c momentum and 
also the tag momentum in the production rate denominator. For example, the upper left entry 

(3.02 ± 0.95) represents the production ratio Ac _[ Ae — I- Ac =2 , for the case where the A c , A c , and 



D each have momentum between 2.3 and 3.3 GeV/c. Statistical errors only are shown. 





2.3 GeV/c< p(Ac) <3.3 GeV/c 


3.3< p(A c ) <5.0 GeV/c 


2.3 GeV/c< p(Kc~ 
3.3 GeV/c< p(A c 


"), p{d" 


) < 3.3 GeV/c 
) <5.0 GeV/c 


3.02 ±0.95 
2.60 ±1.23 


4.04 ± 1.93 
3.79 ± 1.51 


2.3 GeV/c< p(A c ~ 
3.3 GeV/c< p(K c r 


~),p(D~ 


) <3.3 GeV/c 
) <5.0 GeV/c 


3.30 ± 1.33 
3.39 ± 1.76 


3.42 ± 1.99 
3.50 ± 1.47 



The systematic errors are added, mode-by-mode, to determine the overall systematic error for 
each calculated production ratio (Table |v|) . 
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VII. SUMMARY 



The measured values of the production ratios imply that a A c is roughly three times more 
likely to be produced opposite a A c than opposite either a D° or a D~ . These results indicate 
strong evidence in support of correlated production of the A c . Assuming that 6% of charm quark 
fragmentation at ^/s ~10 GeV results in a charmed baryon, our result implies that approximately 
20% of all A c 's are produced in association with a A c in the opposite hemisphere. The effect is 
not produced by the Jetset simulation but does not contradict the Lund String Model. In fact if 
the primary cc were broken at very small proper time by diquark-antidiquark tunneling from the 
vacuum, producing a (possibly excited) charm baryon system and anti-charm antibaryon system, 
a correlation of the type observed here would be generated. 
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